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Abstract 

We present a method developed to deal with electromagnetic wave propagation inside a material 
medium that reacts, in general, non-linearly to the field strength. We work in the context of 
Maxwell's theory in the low frequency limit and obtain a geometrical representation of light paths 
for each case presented. The isotropic case and artificial birefringence caused by an external electric 
field are analyzed as an application of the formalism and the effective geometry associated to the 
wave propagation is exhibited. 

PACS numbers: 04.20.Cv; ll.lO.-z; 42.25.-p 



I. INTRODUCTION 



The characteristics of electromagnetic waves propagating in vacuum are governed by a second order partial differ- 
ential equation. However, when propagating through dispersive media, such differential equation is no longer satisfied 
and, in general, it is necessary to deal with a complicated system of equations. As a consequence, several phenomena 
have to be analyzed closer, like polarization states, wave velocity and its paths inside the particular medium where it is 
traveling. Indeed, such deviation from the linear regime is also found to occur in vacuum when the field strength goes 
beyond its critical value {Ec ^ Be — m^c^/eh), which is the limit of applicability |^] of the linear Maxwell's theory. 
The theoretical description of light propagation in nonlinear electrodynamics was first presented by Bialynicka-Birula 
and Bialynicki-Birula |^ and also by Adler , where the probability of the photon splitting under a strong external 
electromagnetic field was studied. Dittrich and Gies obtained the light cone conditions for a class of homogeneous 
non trivial QED vacua by using the rule of average over polarization states. More recently De Lorenci, Klippert, 
Novello and Salim |^ studied the same problem within a different formalism, in which general light cone conditions 
for a class of theories constructed with the two gauge invariants {F, G) of the Maxwell field was derived without 
making use of average over polarization states. There, birefringence phenomena Q is also described. Further, the 
same authors analyzed the geometrical aspects of wave propagation in the context of nonlinear electrodynamics 
and found that the paths of light can be described in terms of an effective geometry which represents a modification 
of the Minkowski metric. It depends on the dynamics of the background electromagnetic field as 



5r = v"" - 

Lf 



[{Lff + ^±Lfg) F'^^F^" + {Lfg + ^±Lgg) F'^xF^" 



(1) 



where J7± are certain coefficients depending on the dynamics and field strength, and Lx represents the derivative of 
the Lagrangian with respect to the invariant X = (F, G). Indeed, it is shown that photons described by nonlinear 
electrodynamics propagate as null geodesies in such effective metric. (See also Ref. |^ where the effective geometry 
interpretation is considered in the analysis of the Scharnhorst effect.) Therefore, applying the previous formalism 
for material media is a difficult task, since for most cases we do not know an equivalent Lagrangian expression 
describing the physical system. Inside such media, Maxwell's equations are supplemented by the constitutive law 
that relates the electromagnetic excitations D,B and field strengths E,H with the quantities describing the electric 
and magnetic properties of each medium. The constitutive equations are, in general, nonlinear and are determined 
phenomenologically. Schonberg |^ and later Obukhov and Hehl pC| ] derived the equivalent spacetime metric for 
Maxwell's theory on an arbitrary (l-|-3)-dimensional manifold imposing a linear constitutive law. Within the same 
approach, Obukhov, Fukui and Rubilar pT| ] investigated the wave propagation in the Maxwell electrodynamics with 
the most general linear constitutive law and derived the associated Fresnel equation which determines the wave 
normals from the constitutive coefficients. 

In this paper we analyze the electromagnetic wave propagation in material media in the geometrical optics approxi- 
mation. We present the general field equations in terms of some functions depending on the properties of the medium 
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without making any restriction to the constitutive law. The geometrical representation of the light paths for some 
special cases are examined and the effective metric associated to each situation are exhibited. 

In order to achieve more simplicity, in what follows we work in Minkowski spacetime employing a Cartesian coor- 
dinate system. The background metric will be represented by 77^,^, which is defined by diag{+l, — 1, — f , — f ). We set 
the units c = 1. 



II. THE FIELD EQUATIONS 

We define the antisymmetric tensors and P^^ representing the total electromagnetic field. They can be expressed 
in terms of the strengths {E,H) and the excitations {D,B) of the electric and magnetic fields as 

F^, - V^E, - V,E^ - rif^.'^^Vo^Bp (2) 
P^. - V^D, - V,D^ - v,.."^Va,Hp (3) 

where represents the velocity 4-vcctor of an arbitrary observer, which in Galilean coordinates will be given by 
Vf^i = S'^. The Levi-Civita tensor introduced above is defined in such way that r^^i^s _ Since the electric 

and magnetic tensor fields contain only spatial components, we will denote the products as E'^Ea = —E'^ and 
WHa = —H^. In general the properties of the media are determined by the tensors Sa/s and which relate the 
electromagnetic excitation and the field strength by the generalized constitutive laws, 

= £„'3(s^i^^)S0 (4) 

B^ = ^lJ{E^^,H^^)Hp. (5) 

In the absence of sources Maxwell's theory can be summarized by the equations 

V^'D\^~V''D^',^~r|^"'°'^Va,Hp,^^{) (6) 
V^'B^, - VB^^, + r,'""'>^Va,Ep^, = (7) 

which is equivalent to P^'^ = and F^^'^ ^i, = 0, respectively. Therefore, the electromagnetic excitation is related to 
the field strength by means of the constitutive relations (^) and (||), whose derivatives with respect to the coordinates 
can be presented as 

Da,T = £a Ef},r + an FfjEf^^r + o rr EpHf^^r (8) 
Ba.T — Mq Hp T + HpEf^^r + HpH^^r- (9) 

Equations (||) and (0) , together with the relations (^) and , are given in terms of the first derivatives of the electric 
and magnetic fields. They represent the electromagnetic field equations inside an arbitrary material medium whose 
properties are determined by the tensors Sap and /ia/j. 



III. PROPAGATING WAVES 



In order to derive the relations that determine the propagating waves, we will consider the Hadamard method of 
field discontinuities 

iHi- Let us consider a surface of discontinuity E defined by z{x^) — 0. Whenever S is a global 
surface, it divides the spacetime in two distinct regions [/~, for 2: < 0, and f/+, for z > Q. Given an arbitrary function 
of the coordinates, /(x"), we define its discontinuity on S as 



{P±} 



(10) 



where P+, P and P belong to C/+, U 
fields and their derivatives we set 



and S respectively. Applying these conditions for the electric and magnetic 



[E^]^ = 0; [E^_,^y]^ epif,. 



(11) 
(12) 
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where and h^^ represent the discontinuities of the fields on the surface S and K\ is the wave propagation 4-vector. 
Applying these conditions to the field equations (H) and (Q), we obtain the following set of equations governing the 
wave propagation 

e'^^Ko^ep + ——K^Epe^ + ——K^Ef^h,, = (13) 
fi'^^K^hp + -^K^Hpe^ + -^—K^Hphf, = (14) 
V + -Q^Ep^o. + -Q^EphA {KV) + r^^''"'^ K^Vo^hp = (15) 



t^^P + -£^^0^-' + -^HpK {KV) - r^'"'"'' K.Va^ep = (16) 



where we have defined {KV) = K^V^. Indeed, equations (|13| ) and (|lj) follow from the zeroth component of equations 
and (|^), respectively. 

We will restrict our analysis for those cases where Sap = eap{E^ , H'^) and ^afi — ^^'^a^3■ (From the symmetry of 
the field equations, all results we obtain are equally true for the tensor /ia/3). Thus, from equation (14) results the 
identity h^Kp = 0, and equation ( p^ ) gives the following relation between and ea. 

= ^f-'K.V^ep. (17) 

By introducing these results in equation ( p^ ) it follows that 

Z^'^ep = (18) 

with 



Z^'^ = ^i{KV)He>''^ + ^^eA + {KV)^^ijo.^"^EpKrV„ - K^^K^ + {KV)V^Kf^ + if^^^/s _ (ifi/)2^M/3, 
V oEp J oHa 



(19) 



Non-trivial solutions of equation (|l8| ) can be found for such cases where det \ZP'^\ = 0. This condition is known in 
the literature as the Fresnel equation. Despite of its importance, we will not be interested in solving it explicitly, 
but in analyzing some special cases for which an effective metric structure associated to the wave propagation can be 
derived. For the cases where the constitutive laws are linear, the Fresnel equation was solved in Ref. ^0,0. 



IV. THE EFFECTIVE GEOMETRY FOR MEDIA WITH ISOTROPIC PERMITTIVITY 

The special case of media with isotropic permittivity can be stated by assuming e^^ = e{E,H){8^^ — V^V^). 
Introducing this condition in equation (|lq), we obtain 

2 _ (1 _ f,e) {KV)^] 77^'3 „ (i^y)2 1 |i^Mij/3 _ (^y) l_^^r.^0 HrK^V^E^" - Kf^K^ + (if T/)F^X^| ep = 0. (20) 

This equation can be solved by expanding ep as a linear combination of four linearly independent vectors, which can 
be chosen to be Eg, Hp, Kp and as = a^Ep + (igHp + ^oKp + 5aVp. Thus, substituting such representation of 
ep in equation (BQ) and considering the linear independence of the vectors used to expand it, it yields the following 
equations 



K' 



{Kvy 



- 1 + M 



d{tE) _ 1 de 
dE ^ {KV)HdH 



HrK^V^Ep 



EdE 



7o 



^^E-K^ 
EdE 



= 



aoEt^Kp + P.H^Kp + 7o(l - ^ie){KV)^' + 6„{KV) = 

a,{KV)E^Kp + (3„{KV)H^Kp + jo{KV)K^ + So [K^ + ^ie{KV)''] = 

pAK^-{l-lie){KV)''] -0. 



(21) 

(22) 
(23) 
(24) 



The solution of this system results in the light cone condition 
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[KVf 



d{eE) 
dE 



1 dt 

Teoe 



HrK.Va^Ep. 



(25) 



The wave propagation set by equation 
geometry 



5[) can be presented as g'^'^K^Ki, — 0, where we have defined the effective 



III' 
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djeE) 
'^E~ 



V^V^-^^E^E'^-^^ 
eE dE 2H dH 



^-r-^PiPy'^'^HrVa^Ep 



(26) 



in which the notation designates the operation of symmetrization over the index, e.g. 

this way we see that is a null vector in the effective geometry g^^ which constitutes a deviation from Minkowski 
background due to the dielectric media in which the wave propagates. Indeed, it can be shown that /C^, satisfies the 
geodesic equation Kfj_.\K^ = 0, where the semicolon designates the covariant derivative in the geometry determined 
by g^^. Thus, since is a null vector in such effective geometry, the integral curves of will be null geodesies. 

The phase velocity v of the waves can be derived from the light cone condition ( p5| ) by imposing — uP' — 
and defining v = loI\K\, resulting 



djeE) 
' dE 



1 de 
ujH dH 



K-Ex H 



(27) 



where k is an unit vector in the K direction: k = K/\K\. 

In the particular case of homogeneous and isotropic media, where e and fi do not depend on the field strength, the 
effective geometry reduces to the Gordon |13| metric g^"^ = ry'"' + (e/i — l)y^y . The light paths will be determined by 
ds^ — 0, where ds^ = (l/e/i)d<^ — dx^ — rfy^— dz^. It follows that the phase velocity associated to the electromagnetic 
wave is given, as expected, by = 1/ e/i, and the refraction index of the medium Hq — y^e/I. The same result can be 
obtained directly from the general expression (|27|) by imposing e constant. 



V. THE EFFECTIVE GEOMETRY FOR ARTIFICIAL BIREFRINGENCE CAUSED BY AN EXTERNAL 

ELECTRIC FIELD 

When an optically isotropic dielectric is submitted to an external electric field it is experimentally known that it 
may become optically anisotropic - the so called Kerr electro-optic effect. Under such situation the dielectric behaves 
optically like a uniaxial crystal with its optic axis parallel to the field direction . Let us analyze this effect using 
the formalism developed before and derive the effective geometry associated to the wave propagation. 

We set the following relation between the dielectric tensor e^j, and the external field 

£P/3 ^ g „ v'^V^) ~ aE'^E'^ (28) 

with e and a constants. For this case, taking the product of equation (^8|) with E^, and by considering the relations 
( p^ ) and (|l^), we obtain the light cone condition: 

K^=[l-^i{e + 3aE')] {KVf + ^-^^^E^E'^ K^K,. (29) 

Indeed, such condition can be expressed as g^^'''K^Ki, = 0, where we have introduced the effective geometry 

^ +[//(£ + iaE^) - 1] V^V - —^^—r-E^Er (30) 

The discontinuities of the electromagnetic field inside a dielectric that reacts to the action of an external field as 
stated by equation (|2|) , propagates along null geodesies of this effective geometry. 
The phase velocity yields 



2 1 

V — 



^ + ^^(E-kY 
e + aE^ 



(31) 



As before, k is an unit vector in the K direction. We shall recognize two limiting cases from such expression. For 
E ■ k = 0, which leads to vj_ = l/^(e + 3aE^); and for i? • fc = \E\, which leads to v?^ = l/^(e + aE^). 
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In order to identify the birefringence phenomenon we shall examine the eigen-vector problem stated by equation 
( p^ ) . Assuming the initial condition E ■ K = we obtain 

{ [fi{KVy{e + aE^) + ~ (KV)^] 7f ~ 2afi{KVyE^'E''} = 0. (32) 

The solutions for this eigen-vector problem can be achieved by expanding as a linear combination of the linearly 
independent vectors -E^, iJ^, Kf^ and Vp_, resulting in the two light cone conditions on the wave propagation, 

- {KVf [1 - ^l{e + iaE"^)] ^ (33) 
K"^ - {KVf[l- nie + aE"^)] =Q. (34) 

As we see, from equation (p9|), the condition ( ^ ) is already expected. Therefore, we also obtain the second solution 
( ^ ) that corresponds to another mode of light propagation. Performing the computation of the phase velocity 
associated to this mode we obtain ~ 1/M(e + ctE'^), which is the wave existing in the case where E ■ k = \E\ from 
equation (p9[). Such mode survives for all directions of propagation, and there will be an additional mode depending 
on the angle from the wave propagation to the direction of the external electric field. This mode is governed by the 
effective geometry introduced by (|30|). Thus, birefringence occur when an optically isotropic dielectric is submitted 
to an external electric field. 

The effective geometries that govern the electromagnetic wave inside the medium will be given by (^), for the ray 
whose velocity depends on the direction of the external field, and by g^'' = rj'^'^ + [fi{e + aE'^) — l] V^V^, for the ray 
whose velocity does not depend on the direction of the external field. 

For the limiting case, where E ■ K = 0, the phase velocities are determined by and w^, as calculated before and 
the index of refraction yields = /i(e + 3aE^) and = ^{e + aE'^), respectively. In most cases the term in a 
represents only a small correction of the values of e. Thus for electromagnetic waves propagating inside a material 
medium in a perpendicular direction to that defined by the external electric field, the difference between the index of 
refraction of the two existing rays yields 



nj_ — n|| 



aE' a^E^ 



(35) 



where n„ = {efiy^^ is the index of refraction of the isotropic case, as derived in the previous section. Finally when 
a = we recover the isotropic case. 



VI. CONCLUSION 



In this paper we have dealt with electromagnetic wave propagation inside a material medium. The general field 
equations were presented in terms of the field strengths (£^^, i?^) and the tensors e^i, and /x^;^ that characterize the 
properties of each medium where the propagation occurs. By using the Hadamard method we derived the set of 
equations that govern the wave propagation. We examined some special cases where the light cone conditions can 
be obtained without solving the general Fresnel equation and for these cases we presented the associated effective 
geometry. As an application of the formalism, we analyzed the isotropic media and the artificial birefringence caused 
by an external electric field applied over a dielectric medium that reacts nonlinearly to the field excitation. For this 
case we obtained as a first order of approximation the known result that the difference between the index of refraction 
of the two limiting waves propagating in such dielectric is proportional to the square of the external field. 

The method developed here can by applied to study the electromagnetic wave propagation inside material media. 
(We do not considered the frequency dependent material properties.) To do this we must know the expression for the 
tensors e^^ and /i^^ in terms of both the specific medium and the external field. 

An effective geometry can be derived for each situation and can be used in the study of the properties of light 
propagation. With such geometrical description, we present tools for testing kinematic aspects of gravitation in 
laboratory, an issue very much addressed these days (see Refs. p5|-[l^). For instance, we could ask about the 
possibility of formation of structures like event horizons, bending of light and others, in laboratory. 

A possible continuation of this work would be the comparison between the Lagrangian method resumed in the 
effective geometry stated by (|l|) and that obtained for each case that can be derived for non-linear media. Also 
deserves further investigation the general effective geometry associated with an arbitrary material medium in terms 
of the tensors eij,^, and /x^^. 
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